Exploiting microRNA and mRNA profiles generated in vitro from carcinogen-exposed primary mouse hepatocytes for predicting in vivo genotoxicity and carcinogenicity.

The well-defined battery of in vitro systems applied within chemical cancer risk assessment is often characterised by a high false-positive rate, thus repeatedly failing to correctly predict the in vivo genotoxic and carcinogenic properties of test compounds. Toxicogenomics, i.e. mRNA-profiling, has been proven successful in improving the prediction of genotoxicity in vivo and the understanding of underlying mechanisms. Recently, microRNAs have been discovered as post-transcriptional regulators of mRNAs. It is thus hypothesised that using microRNA response-patterns may further improve current prediction methods. This study aimed at predicting genotoxicity and non-genotoxic carcinogenicity in vivo, by comparing microRNA- and mRNA-based profiles, using a frequently applied in vitro liver model and exposing this to a range of well-chosen prototypical carcinogens. Primary mouse hepatocytes (PMH) were treated for 24 and 48h with 21 chemical compounds [genotoxins (GTX) vs. non-genotoxins (NGTX) and non-genotoxic carcinogens (NGTX-C) versus non-carcinogens (NC)]. MicroRNA and mRNA expression changes were analysed by means of Exiqon and Affymetrix microarray-platforms, respectively. Classification was performed by using Prediction Analysis for Microarrays (PAM). Compounds were randomly assigned to training and validation sets (repeated 10 times). Before prediction analysis, pre-selection of microRNAs and mRNAs was performed by using a leave-one-out t-test. No microRNAs could be identified that accurately predicted genotoxicity or non-genotoxic carcinogenicity in vivo. However, mRNAs could be detected which appeared reliable in predicting genotoxicity in vivo after 24h (7 genes) and 48h (2 genes) of exposure (accuracy: 90% and 93%, sensitivity: 65% and 75%, specificity: 100% and 100%). Tributylinoxide and para-Cresidine were misclassified. Also, mRNAs were identified capable of classifying NGTX-C after 24h (5 genes) as well as after 48h (3 genes) of treatment (accuracy: 78% and 88%, sensitivity: 83% and 83%, specificity: 75% and 93%). Wy-14,643, phenobarbital and ampicillin trihydrate were misclassified. We conclude that genotoxicity and non-genotoxic carcinogenicity probably cannot be accurately predicted based on microRNA profiles. Overall, transcript-based prediction analyses appeared to clearly outperform microRNA-based analyses.

Evaluating microRNA profiles reveals discriminative responses following genotoxic or non-genotoxic carcinogen exposure in primary mouse hepatocytes.

Chemical carcinogenesis can be induced by genotoxic (GTX) or non-genotoxic (NGTX) carcinogens. GTX carcinogens have a well-described mode of action. However, the complex mechanisms by which NGTX carcinogens act are less clear and may result in conflicting results between species [e.g. Wy-14,643 (Wy)]. We hypothesise that common microRNA response pathways exist for each class of carcinogenic agents. Therefore, this study compares and integrates mRNA and microRNA expression profiles following short term acute exposure (24 and 48h) to three GTX [aflatoxin B1 (AFB1), benzo[a]pyrene (BaP) and cisplatin (CisPl)] or three NGTX (2,3,7,8-tetrachloordibenzodioxine (TCDD), cyclosporine A (CsA) and Wy) carcinogens in primary mouse hepatocytes. Discriminative gene sets, microRNAs (not for 24h) and processes were identified following 24 and 48h of exposure. From the three discriminative microRNAs found following 48h of exposure, mmu-miR-503-5p revealed to have an interaction with mRNA target gene cyclin D2 (Ccnd2 - 12444) which was involved in the discriminative process of p53 signalling and metabolism. Following exposure to NGTX carcinogens Mmu-miR-503-5p may have an oncogenic function by stimulating Ccnd2 possibly leading to a tumourigenic cell cycle progression. By contrast, after GTX carcinogen exposure it may have a tumour-suppressive function (repressing Ccnd2) leading to cell cycle arrest and to increased DNA repair activities. In addition, compound-specific microRNA-mRNA interactions [mmu-miR-301b-3p-Papss2 (for AFB1), as well as mmu-miR-29b-3p-Col4a2 and mmu-miR-24-3p-Flna (for BaP)] were found to contribute to a better understanding of microRNAs in cell cycle arrest and the impairment of the DNA damage repair, an important hallmark of GTX-induced carcinogenesis. Overall, our results indicate that microRNAs represent yet another relevant intracellular regulatory level in chemical carcinogenesis.

Epigenetic mechanisms underlying arsenic-associated lung carcinogenesis.

Arsenic is an established human carcinogen, but the mechanisms through which it contributes to for instance lung cancer development are still unclear. As arsenic is methylated during its metabolism, it may interfere with the DNA methylation process, and is therefore considered to be an epigenetic carcinogen. In the present study, we hypothesize that arsenic is able to induce DNA methylation changes, which lead to changes in specific gene expression, in pathways associated with lung cancer promotion and progression. A549 human adenocarcinoma lung cells were exposed to a low (0.08 µM), intermediate (0.4 µM) and high (2 µM) concentration of sodium arsenite for 1, 2 and 8 weeks. DNA was isolated for whole-genome DNA methylation analyses using NimbleGen 2.1 M deluxe promoter arrays. In addition, RNA was isolated for whole-genome transcriptomic analysis using Affymetrix microarrays. Arsenic modulated DNA methylation and expression levels of hundreds of genes in a dose-dependent and time-dependent manner. By combining whole-genome DNA methylation and gene expression data with possibly involved transcription factors, a large molecular interaction network was created based on transcription factor-target gene pairs, consisting of 216 genes. A tumor protein p53 (TP53) subnetwork was identified, showing the interactions of TP53 with other genes affected by arsenic. Furthermore, multiple other new genes were discovered showing altered DNA methylation and gene expression. In particular, arsenic modulated genes which function as transcription factor, thereby affecting target genes which are known to play a role in lung cancer promotion and progression.

Elimination of heparin interference during microarray processing of fresh and biobank-archived blood samples.

In the context of environmental health research, biobank blood samples have recently been identified as suitable for high-throughput omics analyses enabling the identification of new biomarkers of exposure and disease. However, blood samples containing the anti-coagulant heparin could complicate transcriptomic analysis because heparin may inhibit RNA polymerase causing inefficient cRNA synthesis and fluorophore labelling. We investigated the inhibitory effect of heparin and the influence of storage conditions (0 or 3 hr bench times, storage at room temperature or -80°C) on fluorophore labelling in heparinized fresh human buffy coat and whole blood biobank samples during the mRNA work-up protocol for microarray analysis. Subsequently, we removed heparin by lithium chloride (LiCl) treatment and performed a quality control analysis of LiCl-treated biobank sample microarrays to prove their suitability for downstream data analysis. Both fresh and biobank samples experienced varying degrees of heparin-induced inhibition of fluorophore labelling, making most samples unusable for microarray analysis. RNA derived from EDTA and citrate blood was not inhibited. No effect of bench time was observed but room temperature storage gave slightly better results. Strong correlations were observed between original blood sample RNA yield and the amount of synthesized cRNA. LiCl treatment restored sample quality to normal standards in both fresh and biobank samples and the previously identified correlations disappeared. Microarrays hybridized with LiCl-treated biobank samples were of excellent quality with no identifiable influence of heparin. We conclude that, to obtain high quality results, in most cases heparin removal is essential in blood-derived RNA samples intended for microarray analysis.

Characterisation of cisplatin-induced transcriptomics responses in primary mouse hepatocytes, HepG2 cells and mouse embryonic stem cells shows conservation of regulating transcription factor networks.

The toxic mechanisms of cisplatin have been frequently studied in many species and in vitro cell models. The Netherlands Toxicogenomics Centre focuses on developing in vitro alternatives using genomics technologies for animal-based assays on, e.g. genotoxic hazards. Models such as human hepatocellular carcinoma cell line (HepG2) cells, mouse primary hepatocytes (PMH) and mouse embryonic stem cells (mESC) are used. Our aim was to identify possibly robust conserved mechanisms between these models using cisplatin as model genotoxic agent. Transcriptomic data newly generated from HepG2 cells and PMH exposed to 7 µM cisplatin for 12, 24 and 48h and 24 and 48h, respectively, were compared with published data from mESC exposed to 5 µM cisplatin for 2-24h. Due to differences in response time between models and marginal changes after shorter exposure periods, we focused on 24 and 48h. At gene level, 44 conserved differentially expressed genes (DEG), involved in processes such as apoptosis, cell cycle, DNA damage response and DNA repair, were found. Functional analysis shows that limited numbers of pathways are conserved. Transcription factor (TF) network analysis indicates 12 common TF networks responding among all models and time points. Four TF, HNF4-α, SP1, c-MYC and p53, capable of regulating ±50% of all DEG, seem of equal importance in all models and exposure periods. Here we showed that transcriptomic responses across several in vitro cell models following exposure to cisplatin are mainly determined by a conserved complex network of 4 TFs. These conserved responses are hypothesised to provide most relevant information for human toxicity prediction and may form the basis for new in vitro alternatives of risk assessment.

Time-series analysis of gene expression profiles induced by nitrosamides and nitrosamines elucidates modes of action underlying their genotoxicity in human colon cells.

N-nitroso compounds (NOCs) may represent a carcinogenic risk to humans following endogenous colonic nitrosation processes. We used the colon adenocarcinoma cell line Caco-2 to investigate transcriptomic changes at three time points (1, 6, 24 h) following exposure to genotoxic concentrations of six different NOCs (two nitrosamides, four nitrosamines) with the purpose of identifying biological processes that may play a part in the carcinogenicity of these compounds. This is especially important for nitrosamide exposure where, in light of their high reactivity, important gene expression modifications may take place early in the exposure. We also analyzed NOC-induced O(6)-methylguanine adducts in relation to transcriptomics since these adducts may influence the expression of genes pivotal in NOC-associated carcinogenicity. Many modified pathways appeared related to DNA damage, cell cycle, apoptosis, growth factor signaling and differentiation, which are linked with carcinogenicity. Nitrosamides showed the strongest response at 1h of exposure, while nitrosamines had the strongest effect at 6 and 24 h. Additionally, methylation was strongly associated with processes that may contribute to the carcinogenic risk. In summary, we have found that NOC-induced gene expression changes vary over time and that many of the modified pathways and processes indicate a carcinogenic risk associated with NOC exposure.

Concentration-dependent gene expression responses to flusilazole in embryonic stem cell differentiation cultures.

The murine embryonic stem cell test (EST) is designed to evaluate developmental toxicity based on compound-induced inhibition of embryonic stem cell (ESC) differentiation into cardiomyocytes. The addition of transcriptomic evaluation within the EST may result in enhanced predictability and improved characterization of the applicability domain, therefore improving usage of the EST for regulatory testing strategies. Transcriptomic analyses assessing factors critical for risk assessment (i.e. dose) are needed to determine the value of transcriptomic evaluation in the EST. Here, using the developmentally toxic compound, flusilazole, we investigated the effect of compound concentration on gene expression regulation and toxicity prediction in ESC differentiation cultures. Cultures were exposed for 24 h to multiple concentrations of flusilazole (0.54-54 µM) and RNA was isolated. In addition, we sampled control cultures 0, 24, and 48 h to evaluate the transcriptomic status of the cultures across differentiation. Transcriptomic profiling identified a higher sensitivity of development-related processes as compared to cell division-related processes in flusilazole-exposed differentiation cultures. Furthermore, the sterol synthesis-related mode of action of flusilazole toxicity was detected. Principal component analysis using gene sets related to normal ESC differentiation was used to describe the dynamics of ESC differentiation, defined as the 'differentiation track'. The concentration-dependent effects on development were reflected in the significance of deviation of flusilazole-exposed cultures from this transcriptomic-based differentiation track. Thus, the detection of developmental toxicity in EST using transcriptomics was shown to be compound concentration-dependent. This study provides further insight into the possible application of transcriptomics in the EST as an improved alternative model system for developmental toxicity testing.

Evaluation of developmental toxicant identification using gene expression profiling in embryonic stem cell differentiation cultures.

The murine embryonic stem cell test (EST) is an alternative testing method designed to assess potential developmental toxicity of compounds. The implementation of transcriptomics in the EST has been shown to reduce the culture duration and improve endpoint evaluation and is expected to result in an enhanced predictability and definition of the applicability domain. We evaluated the identification of developmental toxicity in the EST using two gene sets ("Van_Dartel_heartdiff_24h" and "EST biomarker genes") defined in our earlier studies. Nonexposed embryonic stem cells (ESC) differentiation cultures were sampled 0, 24, and 48 h after initiation of differentiation. Additionally, cultures exposed to 12 diverse well-characterized positive and negative developmental toxicants were isolated 24 h after the onset of exposure. Inhibition of ESC differentiation was evaluated in parallel by morphological scoring on culture day 10. Transcriptomics analysis was conducted using the Affymetrix Gene Chips platform. We applied principal component analysis on the basis of the two predefined gene sets to define the "differentiation track" that represents ESC differentiation. The significance of derivations in the gene expression-based differentiation track because of compound exposures were evaluated to determine developmental toxicity of tested compounds. We successfully predicted developmental toxicity using transcriptomics for 83% (10/12) and 67% (8/12) of the compounds, respectively, using the two predefined gene sets ("Van_Dartel_heartdiff_24h" and "EST biomarker genes"). Our study suggests that the application of transcriptomics may improve the applicability of the EST for the prediction of the developmental toxicity of chemicals.

Time series analysis of benzo[A]pyrene-induced transcriptome changes suggests that a network of transcription factors regulates the effects on functional gene sets.

Chemical carcinogens may cause a multitude of effects inside cells, thereby affecting transcript levels of genes by direct activation of transcription factors (TF) or indirectly through the formation of DNA damage. As the temporal profiles of these responses may be profoundly different, examining time-dependent changes may provide new insights in TF networks related to cellular responses to chemical carcinogens. Therefore, we investigated in human hepatoma cells gene expression changes caused by benzo[a]pyrene at 12 time points after exposure, in relation to DNA adduct and cell cycle. Temporal profiles for functional gene sets demonstrate both early and late effects in up- and downregulation of relevant gene sets involved in cell cycle, apoptosis, DNA repair, and metabolism of amino acids and lipids. Many significant transcription regulation networks appeared to be around TF that are proto-oncogenes or tumor suppressor genes. The time series analysis tool Short Time-series Expression Miner (STEM) was used to identify time-dependent correlation of pathways, gene sets, TF networks, and biological parameters. Most correlations are with DNA adduct levels, which is an early response, and less with the later responses on G1 and S phase cells. The majority of the modulated genes in the Reactome pathways can be regulated by several of these TF, e.g., 73% by nuclear factor-kappa B and 34-42% by c-MYC, SRF, AP1, and E2F1. All these TF can also regulate one or more of the others. Our data indicate that a complex network of a few TF is responsible for the majority of the transcriptional changes induced by BaP. This network hardly changes over time, despite that the transcriptional profiles clearly alter, suggesting that also other regulatory mechanisms are involved.

Discrimination for genotoxic and nongenotoxic carcinogens by gene expression profiling in primary mouse hepatocytes improves with exposure time.

Assessing the potential carcinogenicity of chemicals for humans represents an ongoing challenge. Chronic rodent bioassays predict human cancer risk at only limited reliability and are simultaneously expensive and long lasting. In order to seek for alternatives, the ability of a transcriptomics-based primary mouse hepatocyte model to classify carcinogens by their modes of action was evaluated. As it is obvious that exposure will induce a cascade of gene expression modifications, in particular, the influence of exposure time in vitro on discriminating genotoxic (GTX) carcinogens from nongenotoxic (NGTX) carcinogens class discrimination was investigated. Primary mouse hepatocytes from male C57Bl6 mice were treated for 12, 24, 36, and 48 h with two GTX and two NGTX carcinogens. For validation, two additional GTX compounds were studied at 24 and 48 h. Immunostaining of gammaH2AX foci was applied in order to phenotypically verify DNA damage. It confirmed significant induction of DNA damage after treatment with GTX compounds but not with NGTX compounds. Whole-genome gene expression modifications were analyzed by means of Affymetrix microarrays. When using differentially expressed genes from data sets normalized by Robust Multi-array Average, the two classes and various compounds were better separated from each other by hierarchical clustering when increasing the treatment period. Discrimination of GTX and NGTX carcinogens by Prediction Analysis of Microarray improved with time and resulted in correct classification of the validation compounds. The present study shows that gene expression profiling in primary mouse hepatocytes is promising for discriminating GTX from NGTX compounds and that this discrimination improves with increasing treatment period.

Assessing the metabolic competence of sandwich-cultured mouse primary hepatocytes.

Primary human and rat hepatocyte cultures are well established in vitro systems used in toxicological studies. However, whereas transgenic mouse models provide an opportunity for studying mechanisms of toxicity, mouse primary hepatocyte cultures are less well described. The potential usefulness of a mouse hepatocyte-based in vitro model was assessed in this study by investigating time-dependent competence for xenobiotic metabolism and gene expression profiles. Primary mouse hepatocytes, isolated using two-step collagenase perfusion, were cultured in a collagen sandwich configuration. Gene expression profiles and the activities of various cytochrome P450 (P450) enzymes were determined after 0, 42, and 90 h in culture. Principal component analysis of gene expression profiles shows that replicates per time point are similar. Gene expression levels of most phase I biotransformation enzymes decrease to approximately 69 and 57% of the original levels at 42 and 90 h, respectively, whereas enzyme activities for most of the studied P450s decrease to 59 and 34%. The decrease for phase II gene expression is only to 96 and 92% of the original levels at 42 and 90 h, respectively. Pathway analysis reveals initial effects at the level of proteins, external signaling pathways, and energy production. Later effects are observed for transcription, translation, membranes, and cell cycle-related gene sets. These results indicate that the sandwich-cultured primary mouse hepatocyte system is robust and seems to maintain its metabolic competence better than that of the rat hepatocyte system.

Development and validation of a modified comet assay to phenotypically assess nucleotide excision repair.

There is an increasing need for simple and reliable approaches to phenotypically assess DNA repair capacities. Therefore, a modification of the alkaline comet assay was developed to determine the ability of human lymphocyte extracts to perform the initial steps of the nucleotide excision repair (NER) process, i.e. damage recognition and incision. Gel-embedded nucleoids from A549 cells, pre-exposed to 1 microM benzo[a]pyrene-diol-epoxide, were incubated with cell extracts from frozen or freshly isolated lymphocytes. The rate at which incisions are introduced and the subsequent increase in tail moment is indicative for the repair capacity of the extracts. Freshly prepared extracts from lymphocytes of human volunteers (n = 8) showed significant inter-individual variations in their DNA repair capacity, which correlated with the removal of bulky DNA lesions over a period of 48 h determined by (32)P-post-labelling (R(2) = 0.76, P = 0.005). Repeated measurements revealed a low inter-assay variation (11%). Storage of cell extracts for more than 3 weeks significantly reduced (up to 80%) the capacity to incise the damaged DNA as compared to freshly isolated extracts. This reduction was completely restored by addition of ATP to the extracts before use, as it is required for the incision step of NER. In contrast, extracts freshly prepared from frozen lymphocyte pellets can be used without loss of repair activity. DNA repair deficient XPA-/- and XPC-/- fibroblasts were used to further validate the assay. Although some residual capacity to incise the DNA was observed in these cells, the repair activity was restored to normal wild-type levels when a complementary mixture of both extracts (thereby restoring XPA and XPC deficiency) was used. These results demonstrate that this repair assay can be applied in molecular epidemiological studies to assess inter-individual differences in NER.